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1.  INTRODUCTION 


In  order  to  achieve  accurate  artillery  fire,  the  net  effect  of  winds  on  the 
projectile  at  various  altitudes  must  be  taken  into  consideration.  One  aspect 
of  quantifying  this  effect  involves  the  computation  of  ballistic  wind  velocity 
from  data  collected  by  any  of  several  types  of  meteorological  (MET)  data 
acquisition  systems.  Inherent  errors  in  measurement  propagate  to  produce  an 
undesirable  error  in  the  computed  value  of  ballistic  wind,  since  the  measure¬ 
ment  errors  depend  on  the  particular  type  of  MET  system  employed.  It  is 
useful  to  have  methods  by  which  the  performance  of  various  systems  of  interest 
can  be  compared. 

In  this  report,  three  error  analysis  models  are  presented,  each  of  which  is 
applicable  to  one  general  type  of  MET  data  system.  Each  system  modeled  mea¬ 
sures  a  particular  set  of  variables  to  determine  successive  positions  of  a 
balloon-borne  radiosonde  as  it  ascends  through  the  atmosphere.  Known  or 
estimated  errors  in  the  geometric  variables,  an  actual  or  postulated  wind 
profile,  and  certain  other  parameters  are  inputted  to  the  appropriate  model. 
From  this  information  the  model  computes,  among  other  results,  a  quantity 
called  the  component  velocity  variance  (CVV)  of  the  ballistic  wind. 

For  a  given  direction  of  artillery  fire,  the  ballistic  wind  can  be  resolved 
into  range  wind  and  crosswind  components.  The  CVV  is  obtained  by  averaging 
the  variance  associated  with  either  of  these  components  over  all  possible 
directions  of  fire  in  the  horizontal  plane.  The  CVV  represents  a  general 

error  quantity  that  can  be  used  as  a  basis  for  comparing  MET  systems  or  for 
evaluating  suggested  Improvements  in  them. 

The  following  general  assumptions  are  made:  (I)  the  CVV  exists  and  Is  a 

useful  characterization  of  the  measuring  system,  and  (2)  first-order  error 
analysis  is  sufficient  to  determine  the  CVV  to  an  acceptable  accuracy. 

Although  not  absolutely  necessary,  it  Is  convenient  to  assume  that  errors  In 

measurement  are  normally  distributed.  Further  assumptions  are  noted  later  as 
they  are  required. 


The  error  analysis  models  are  outlined  below: 

(1)  RAW1N  models  the  case  of  a  balloon-borne  radiosonde  and  a  ground-based 
set  for  radiodirection  finding  (RDF)  and  telemetry  data  reception.  This 
model  is  applicable,  for  example,  to  the  Rawin  Set  AN/GMD-1.  It  is  also 
applicable  to  the  Meteorological  Data  System  AN/TMQ-31,  operating  in  the 
RDF  mode;  this  set  is  also  known  as  the  Field  Artillery  Meteorological 
Acquisition  System  (FAMAS). 

(2)  RADAR  models  the  case  in  which  a  ground-based  radar  measures  all  the 
variables  required  to  determine  the  balloon's  position.  This  model  is 
applicable,  for  example,  to  the  radar  sets  AN/TMQ-19  and  AN/FPS-16. 

(3)  NAVAID  models  the  case  in  which  radionavigation  techniques  are  used  to 
determine  the  position  of  the  radiosonde.  This  model  is  applicable  to 
the  radionavigation  portions  (LORAN,  VLF,  and  OMEGA)  of  the  AN/TMQ-31. 

The  models  represent  a  compromise  between  extreme  generality  and  overt  spe¬ 
cialization  to  any  particular  field  environment. 

Material  relevant  to  all  three  models  is  presented  in  Sections  2,  3,  4,  and  fi 
of  this  report.  Sections  5,  6,  and  7  discuss  specialization  to  each  of  the 
models,  respectively.  Part  of  the  mathematical  treatment  Is  relegated  to  the 
appendices.  Computer  programs  and  sample  calculations  are  given  in  the  accom¬ 
panying  User's  Manual. 
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2.  PRELIMINARY  CONSIDERATIONS 


For  the  purpose  of  computing  ballistic  wind,  the  atmosphere  is  divided  into  a 
series  of  zones,  not  all  of  which  are  necessarily  of  the  same  vertical  thick¬ 
ness.  Ballistic  line  n  is  associated  with  the  path  followed  by  a  round  that 
attains  its  maximum  altitude  at  the  top  of  the  n'th  zone.  The  components  of 
ballistic  wind  for  line  n  can  be  determined  from  the  average  wind  components 
in  each  of  the  first  n  zones  in  conjunction  with  a  set  of  zone  wind  weighting 
factors  appropriate  to  line  n.  Background  information  on  the  various  accepted 
zone  structures  and  weighting  factors  is  available  in  References  1  and  2. 

In  either  firing  table  or  computer  gunnery,  the  weighting  factors  given  in 
Reference  2  are  believed  to  be  a  good  representation  of  meteorological  effects 
on  commonly  used  surface-to-surface  and  surface-to-air  ballistic  projec¬ 
tiles.  However,  should  the  user  desire,  he  may  readily  substitute  other 
appropriate  weighting  factors,  which  can  be  computed  from  results  of  ballistic 
simulation  programs. 

To  obtain  the  data  required  to  compute  the  ballistic  wind,  a  balloon-borne 
radiosonde  is  tracked.  In  the  error  models  presented  here,  it  is  assumed  that 
the  radiosonde,  actual  or  postulated,  ascends  at  a  constant  rate.  It  is  also 
assumed  that  the  average  velocity  of  the  wind  in  any  zone  is  the  same  as  the 
average  horizontal  velocity  of  the  radiosonde  package  as  it  traverses  the 
zone.  This  is  reasonable  because  of  the  relatively  large  area  of  a  MET  bal¬ 
loon  and  the  relatively  small  mass  of  a  balloon-radiosonde  combination.  As  is 
customary  among  physicists,  engineers,  and  mathematicians,  we  take  the  wind 
direction  to  be  that  toward  which  the  wind  is  blowing;  this  sense  is  opposite 
to  the  usual  meteorological  convention.  (The  computed  values  of  the  CVV  are 
not  affected  by  the  choice  of  wind  convention.) 

Figure  l  illustrates  the  geometric  variables  associated  with  the  radiosonde's 
position  at  the  top  of  the  i'th  ballistic  zone.  In  the  RAWIN  and  RADAR 
models,  the  ground-based  RAWIN  or  RADAR  set  defies  the  origin  of  coordinates; 
in  NAVAID  an  arbitrary  origin  can  be  chosen.  At  time  tj  the  radiosonde  is  at 
altitude  z^  at  the  top  of  zone  i;  its  elevation  angle  is  and  its  azimuth 
angle  is  ajj  and  y ^  are,  respectively,  the  East  and  North  components  of  the 
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Figure  1.  Geometric  Variables 


East 


radiosonde's  position;  is  the  distance  along  the  surface  of  the  earth  from 
the  origin  to  a  point  directly  below  the  radiosonde;  Sj  is  the  slant  range. 
For  ballistic  line  n,  the  zone  index  1  takes  on  the  values  1,  2,  .  ..,  n;  the 
value  i  =  0  corresponds  to  the  launch  position  and  is  treated  as  a  special 
case.  The  set  of  variables  taken  as  independent  (or  even  utilized)  depends  on 
the  particular  model.  For  simplicity  of  presentation.  Figure  1  illustrates 
the  variables  for  the  case  of  a  flat  earth;  however,  in  all  three  error  analy¬ 
sis  models,  computations  are  carried  out  assuming  a  spherical  earth. 


All  three  models  calculate  the  variance  associated  with  each  component  (East 
and  North)  of  ballistic  wind  prior  to  finding  the  CVV.  These  variances  are 
computed  from  the  following  inputs: 


(1)  The  zone  structure  and  appropriate  weighting  factors 

(2)  The  ascent  rate  of  the  balloon 

(3)  The  zone  wind  profile 

(4)  The  successive  values  of  the  Independent  geometric  variable  determining 
the  position  of  the  radiosonde 

(5)  The  errors  associated  with  these  variables 

The  input  errors  are  taken  to  be  standard  deviations.  Explanations  of  all 
required  inputs  are  tabulated  in  the  sections  of  this  report  that  deal  with 
the  individual  models. 

In  RAWIN  and  NAVAID,  errors  associated  with  the  determination  of  altitude  must 
be  estimated  beforehand  or  calculated  from  actual  pressure  and  temperature 
data;  no  provision  exists  in  the  models  for  dealing  directly  with  pressure  and 
temperature  measurements.  See  Reference  3,  for  example. 

In  RAWIN  and  RADAR,  errors  associated  with  the  measurement  of  elevation  angle 
are  divided  into  two  groups:  (1)  errors  In  elevation  due  to  reflection  of  the 
incoming  signal  at  the  surface  of  the  earth,  and  (2)  errors  in  elevation  due 
to  all  other  causes.  The  latter  group  Is  referred  to  as  the  errors  in  the 
elevation  tracking  of  the  apparent  target;  for  brevity  it  is  also  referred  to 
here  simply  as  "leveling"  errors.  The  leveling  errors  can  be  estimated  or 
obtained  from  equipment  manufacturer's  specifications,  while  the  reflection 
errors  must  be  estimated  or  computed  externally  to  the  models  for  given  soil 
types  and  signal  frequencies. 

It  is  assumed  that  any  uncertainty  in  the  actual  measurement  .  time  is  negli¬ 
gible.  However,  because  t^  represents  the  elapsed  time  as  the  balloon  passes 
the  top  of  zone  i,  it  contains  some  error  because  there  is  error  in  the  deter¬ 
mination  of  altitude.  Since  the  balloon  has  a  constant  ascent  rate,  the  error 
in  t^  is  directly  related  to  the  error  in  altitude.  In  actuality,  the  models 
use  time  only  implicitly  as  a  parameter  that  establishes  a  correspondence 
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among  the  geometric  variables.  In  RAWIN,  for  example,  this  simply  means  that 
values  of  altitude,  elevation,  and  azimuth  are  available  as  the  radiosonde 
passes  each  zone  top. 

Two  general  categories  of  error  are  treated:  (1)  bias  errors,  which  are 
correlated  from  zone  to  zone;  and  (2)  random  errors,  which  are  uncorrelated 
from  zone  to  zone.  These  are  discussed  further  in  Section  4. 

The  models  are  programmed  to  accept  the  NATO  zone  structure  of  15  ballistic 
zones  as  input,  along  with  the  necessary  zone  wind  weighting  factors. 
However,  the  models  can  be  modified  to  accept  an  arbitrary  zone  structure 
containing  up  to  30  zones,  along  with  appropriate  weighting  factors.  This 
modification  is  explained  in  the  accompanying  Users'  Manual.  Also,  each  model 
will  accept  up  to  four  balloon  ascent  rates  simultaneously. 

In  performing  Intermediate  computations,  each  model  converts  all  inputs  to 
meter-kilogram-second  (MRS)  units.  Many  of  the  Intermediate  results  are 
optionally  available  as  output.  The  CVV  itself  is  computed  for  each  ballistic 
line  for  each  balloon  ascent  rate  and  is  outputted  in  knots  squared.  The 
square  root  of  the  CVV  is  called  the  standard  deviation  and  is  outputted  in 
knots. 

In  the  following  sections  of  this  report,  various  equations  required  by  the 
models  are  derived  for  the  case  of  a  single  ballistic  line  and  a  single 
balloon  ascent  rate.  Unless  otherwise  noted,  all  variables  and  constants  used 
in  these  derivations  are  considered  to  be  in  MRS  units.  Extension  to 
additional  ballistic  lines  and  ascent  rates  is  straightforward. 
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3.  THE  BALLISTIC  WIND 


The  three  error  models  described  in  thiB  report  make  common  use  of  certain 
algebraic  expressions  Involving  the  components  of  ballistic  wind.  In  this 
section,  these  expressions  are  derived  for  the  East  component.  Equations 
applicable  to  the  North  component  can  be  obtained  in  an  analogous  fashion.  In 
all  the  equations,  it  is  implicitly  recognized  that  the  zone  index  i  takes  on 
all  values,  i  ■  1,  2,  ....  n,  appropriate  to  ballistic  line  n.  The  launch 
index,  i.e.,  i  *  0,  is  treated  separately  as  a  special  case. 

Let  uj  and  v^  be,  respectively,  the  East  and  North  components  of  the  average 
wind  velocity  in  zone  1.  They  are  defined  by  the  expressions 


and 


(3-1) 


(3-2) 


where  x^  and  y^  are,  respectively,  the  East  and  North  coordinates  of  the 
radiosonde  at  time  t^  at  the  top  of  zone  i. 

The  ballistic  wind  is  that  single  value  that  is  equivalent  to  the  cumulative 
effect  of  the  individual  zone  winds.  The  East  component  U  of  the  ballistic 
wind  is  defined  by: 
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where  Wj-  is  the  wind  weighting  factor  appropriate  to  the  i'th  zone  for  ballis¬ 
tic  line  n.  Then, 
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where  Xq  is  the  East  coordinate  of  launch  position,  and  t^  is  the  launch  time. 
Collecting  terms  in  yields 
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where  w  ..  =  0  for  ballistic  line  n. 
n+1 

For  a  constant  ascent  rate  v  ,  the  time  spent  by  the  balloon  in  zone  i  is 
given  by: 
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where  z^  is  the  altitude  at  the  top  of  the  i'th  zone.  Combining  Eq.  (3-7) 
with  Eq.  (3-6),  we  have 


U 


v 

z 


Z1  -  Z0 


X0  + 


a 

I 

i= 1 


w . 
l 


2 

1 


2 


i-1 


w  i+i 


'i+l 


—  1  x. 

Z  .  1 


(3-8) 


14 


The  set  of  altitudes  z^  defines  the  2one  structure  and  is  available  in  the 
form  of  tabulated  values.  From  these  values,  it  is  convenient  to  define  for 
ballistic  line  n  a  set  of  new  constant  weighting  factors  (per  unit  length),  Wq 
and  W^,  that  depend  only  on  the  zone  structure: 


Z1  ‘  20 


(3-9) 
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Then  Eq.  (3-8)  becomes 
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(3-11) 


8U 

The  partial  derivatives  will  be  required  to  compute  the  variance  in  the  East 

i 

component  of  ballistic  wind.  From  Eq.  (3-11)  these  are: 


au 

5 - =  v  Wn 

3Xq  z  0 


(3-12). 


and 


au 

5 —  =  V  W. 

3x .  z  l 

l 


(3-13) 


au 

The  derivatives  are  also  required.  To  compute  these,  we  consider  only 

those  terms  of  Eq .  (3-4)  that  contain  the  general  zone  index  i  explicitly: 


w. (x.  -  X .  . )  w. , . (x.  -  X. ) 

U  =  ...  +  -i— j— — 111!  +  - - L_  + 


t.  -  t.  . 
l  i-l 


t ...  -  t. 
l  +  l  l 


(3-14) 
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Then 


3U 

at. 

i 


wi(xi  ~  xi-i)  .  wi+i(xi+i  *  xi} 
(ti  -  Vi^  (ti+i  -  V*' 


(3-15) 


Making  use  of  the  relationships  expressed  in  Eqs .  (3-1)  and  (3-7),  we  rewrite 
Eq.  (3-15)  as: 


3 U  _  ,  Wi+1 

3t .  z.,  -  z.  ui  +  l 

1  1+1  1 


w . 
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z .  -  z .  _ 

1  1-1 


u .  ]  v 

1  z 


(3-16) 


31) 

We  could  also  find  5 — ,  but  an  explicit  determination  of  this  quantity  is  not 
3t0 

required  by  the  error  analysis  models. 

From  Eq.  (3-7)  we  have 


3t. 

3z .  Vz 

1 


(3-17) 


3U‘  9t. 

The  models  make  use  of  the  product  of  derivatives 

and  (3-17)  we  have  1  1 


From  Eqs.  (3-16) 
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1  1 


Vi 


"i+1 


z .  i+1 
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u . 
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(3-18) 


For  convenience  we  write  this  more  briefly  as 


31)  3t 

3t .  3z . 
1  1 


~  W 


ui 


(3-19) 


where 
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(3-20) 
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Similar  expressions  involving  V,  the  North  component  of  ballistic  wind,  can 
also  be  derived.  The  models  make  use  of  the  following: 


3V  u 

3y0  vz  0  ’ 


3V  u 

~ —  =  v  W.  , 

3y.  z  1 


and 


av 


at. 


at.  3z. 


w  .  , 
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w 

i+1 


v . 


w . 
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vi  z . . ,  -  z .  i+1  z.-z..  i 

l+l  l  l  l-l 


(3-21) 


(3-22) 


(3-23) 


(3-24) 
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4.  FORM  OF  THE  CVV 


The  variance  in  either  component  of  ballistic  wind  can  be  written  in  a  general 
form  that  is  easily  specialized  to  each  of  the  error  models  under  considera¬ 
tion.  In  this  section,  we  examine  this  form  and  the  associated  assumptions 
concerning  bias  and  random  error.  Finally,  we  obtain  the  expression  for  the 
CVV. 


From  Eq.  (3-4),  the  East  component  U  of  ballistic  wind  for  ballistic  line  n 
can  be  written  in  the  implicit  functional  form 


U  =  U  {xQ,  tQ,  x.,  tiJ.  i  =  1,  2,  ...  n  . 

In  general,  each  of  the  variables  Xq,  Cq,  x^,  and  t.  may  themselves  be  func¬ 
tions  of  other  variables.  In  RAWIN,  for  example,  the  East  component  Xq  of 
launch  position  is  taken  to  be  a  function  of  Do  and  cxq,  where  Do  is  the 
distance  from  the  Rawin  Set  to  the  launch  site,  and  (Xq  is  the  launch  site  azi¬ 
muth.  Also  in  RAWIN,  the  ascending  radiosonde's  East  component  x^  is  consid¬ 
ered  to  be  a  function  of  the  altitude  z^,  elevation  Cj,  and  azimuth  a^,  all 
determined  for  the  top  of  zone  i.  In  all  the  models,  tj  is  taken  to  be  a 
functioi  of  Zj  through  Eq.  (3-7). 

In  general,  we  can  rewrite  Eq .  (4-1)  as 

U  =  1)  f  A  j  q  ,  ^20’  \o  ’  ^li  ’  ^2i  ’  -  n  >  (^”2) 

where  X^q,  ^20*  ^L0  the  ra°del-dependent  subset  of  independent  launch 

variables,  and  •  ••,  £]<£  comprise  the  model-dependent  subset  of  inde¬ 

pendent  variables  appropriate  to  zone  i. 

For  a  given  independent  variable,  two  categories  of  measurement  error  are 
considered:  random  error  and  bias  error.  The  random  errors  of  measurement 

follow  some  distribution  law  (which  for  convenience  may  be  taken  to  be  normal) 
with  zero  mean  and  a  characteristic  variance  and  standard  deviation. 


The  bias  error  in  a  given  variable  may  arise  from  inaccurate  calibration  of 
the  measuring  device,  causing  the  determined  values  to  be  always  too  low  or 
too  high  for  a  particular  experiment.  However,  over  many  experiments,  the 
calibration  of  an  instrument  is  as  likely  to  be  overper formed  as  under- 
performed.  From  this  standpoint,  the  bias  errors  themselves  are  random  and 
have  some  distribution  (which  for  convenience  may  be  assumed  to  be  normal) 
with  zero  mean  and  a  characteristic  variance  and  standard  deviation.  This  is 
the  point  of  view  adopted  in  this  report.  Bias  and  random  errors  are  treated 
differently  here  only  in  the  assumptions  made  concerning  their  zone-to-zone 
correlations. 

We  expect  no  correlation  among  the  variables  4^1  and  4ji>  1  “  1.  2,  ....  n, 
j  /  k,  because  values  of  these  variables  are  determined  from  completely  dif¬ 
ferent  types  of  measurement.  For  example,  4li>  42i>  anc*  ^3i  m*8ht  represent 
zl>  el>  and  ai>  respectively. 

It  is  assumed  that  errors  in  the  measurement  of  any  one  type  of  zone  variable 
£ki»  i  =  1,  2,  ...,  n  can  be  treated  in  the  following  manner:  random  errors 
are  uncorrelated  from  zone  to  zone;  bias  errors  are  completely  correlated  from 
zone  to  zone.  In  other  words,  the  random  errors  have  correlation  coefficients 
of  zero,  while  the  bias  errors  have  correlation  coefficients  of  one. 

Finally,  we  assume  only  random  errors  in  the  determination  of  the  independent 
launch  variables  Aio>  ^20«  •••* 


If  only  first-order  propagation  of  error  Is  considered,  it  is  shown  in  Appen¬ 
dix  A  that  the  variance  cru2  la  the  East  component  of  ballistic  wind  for  line  n 
is  given  by: 


au 

aXTT  °RAkO 
kO 


2  K 
]  +  I 
k=l 


2  K 
]  +  I 
k=l 


n 

I 

i=l 


(4-3) 
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In  Eq.  (4-3),  is  the  standard  deviation  associated  with  random  errors  in 
the  measurement  of  Ako.  while  OB^ki  and  OR£ki  are  the  standard  deviations 
associated  respectively  with  the  bias  and  random  errors  in  the  measurement  of 
£ki*  In  computation,  the  bias  contributions  are  summed  over  the  appropriate 
zones  before  being  squared,  while  the  random  contributions  are  squared  before 
being  summed.  Similarly,  the  variance  ay2  in  the  North  component  V  of  ballis¬ 
tic  wind  is 
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k=l 


3V 

3\,  .  °RAkO 
kO 


2  K 

]  +  I 

k=l 


n 

I 

i=l 


3V  „ 

a4ki  B^ki 


2  K 
]  ♦  2 
k=l 
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1 

i=l 


3V 


2 
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(4-4) 


For  a  given  direction  of  artillery  fire,  the  ballistic  wind  can  be  resolved 
into  range  wind  and  crosswind  components.  The  CVV  is  obtained  by  averaging 
the  variance  associated  with  either  of  these  components  over  all  possible 
directions  of  fire  in  the  horizontal  plane. 

Let  the  j'th  direction  of  artillery  fire  make  an  angle  8j  with  the  positive  X 
axis  or  East  direction.  We  define  in  the  horizontal  plane  an  orthogonal 
coordinate  system,  labeled  XjYj,  such  that  the  positive  Xj  axis  points  in  the 
(horizontal)  direction  of  fire.  This  coordinate  system  is  depicted  in  Figure 
2.  Each  possible  direction  of  fire  corresponds  to  a  particular  orientation  of 
the  XjYj  system. 

The  Xj  component  of  ballistic  wind  is  the  range  wind  Uj,  and  the  Yj  component 
is  the  crosswind  Vj,  as  shown  in  Figure  2.  From  a  standard  transformation  of 
coordinates,  it  can  be  shown  that 


U.  =  U  cos  0.  +  V  sin  0. 
J  J  J 


(4-5) 
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(4-12) 


Each  of  the  first  two  integrals  yields  7t,  while  the  third  integral  yields 
zero.  Then 


(4-13) 


The  quantity  oc2  Is  the  CVV  that  is  referred  to  elsewhere  In  this  report.  It 
is  used  to  characterize  the  error  behavior  of  each  MET  data  system. 
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5.  RAWIN 


RAWIN  models  the  case  of  a  balloon-borne  radiosonde  and  a  ground-based  set  for 
RDF  and  telemetry  data  reception.  The  angular  variables  of  the  radiosonde's 
position  are  measured  by  the  RDF  ground  set,  and  the  altitude  can  be 
determined  from  MET  messages  received  from  the  radiosonde.  This  model  is 
applicable,  for  example,  to  the  RAWIN  set  AN/GMD-1  and  to  the  MET  data  system 
AN/TMQ-31  (FAMAS)  operating  in  the  RDF  mode. 

In  RAWIN,  the  radiosonde's  East  coordinate  is  assumed  to  be  a  function  of 
the  independent  variables  of  altitude  z^,  elevation  angle  e^,  and  azimuth 
angle  a^,  all  determined  for  the  top  of  the  i'th  ballistic  zone.  These 
variables  are  illustrated  in  Figure  1.  The  time  t^  at  which  the  radiosonde 
passes  the  top  of  the  zone  is  taken  to  be  a  function  of  zj.  The  balloon  is 
released  at  time  tg.  The  East  coordinate  xq  of  launch  position  is  taken  to  be 
a  function  of  Dq  and  aQ,  where  Dq  is  the  distance  from  the  receiving  set  to 
the  launch  site,  and  oq  is  the  associated  azimuth. 

Under  the  conditions  noted  above,  Eq .  (4-1)  may  be  written  in  the  slightly 

more  explicit  functional  form: 

U  =  U  [k0  (D0,  «„)'  t0,  x.  (z,,  ct,  (r),  tj  (2,)!.  i  =  1,  2,  ... 

(5-1) 

where  U  is  the  East  component  of  ballistic  wind  appropriate  to  line  n.  A 
similar  equation  holds  for  the  North  component. 

To  obtain  the  variance  associated  with  U,  we  rewrite  Eq .  (4-3)  in  terms  of 
the  independent  variables  Dq,  Oq,  z^,  e^,  and  a^.  This  yields 
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In  Eq .  (5-2),  each  type  of  error  is  designated  by  o  with  appropriate  sub¬ 
scripts.  The  first  subscript,  B  or  R,  identifies  the  error  as  bias  or  random, 
respectively.  The  second  subscript  is  the  variable  containing  the  error.  The 
third  subscript  is  0  for  the  launch  variables;  otherwise,  it  is  the  zone  index 
i.  For  errors  in  the  measurement  of  elevation,  a  fourth  subscript  is  pres¬ 
ent:  G  identifies  the  error  associated  with  the  reflection  of  the  transmitted 
signal  by  the  ground;  L  identifies  the  combined  error  in  elevation  due  to  all 
other  causes.  Thus,  for  example,  is  the  random  error  in  the  measurement 
of  elevation  at  zone  top  i  due  to  ground  reflection.  It  is  assumed  that  the 
launch  time  is  known  precisely;  hence,  an  error  term  involving  tg  is  not 
included  in  Eq .  (5-2). 

It  is  assumed  that  elevation  errors  due  to  ground  reflection  are  not  corre¬ 
lated  with  elevation  errors  due  to  other  causes.  Thus  they  enter  indepen¬ 
dently  into  the  computation  of  og2  in  the  manner  shown  in  Eq .  (5-2). 

Analogous  to  Eq.  (5-2),  the  expression  for  the  variance  Oy2  in  the  North 
component  V  of  ballistic  wind  is 
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(5-3) 
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In  order  to  obtain  computationally  useful  error  coefficients,  each  of  the 
partial  derivatives  in  Eqs.  (5-2)  and  (5-3)  must  be  further  expanded.  The 
form  of  Eq.  (5-1)  allows  for  straightforward  use  of  the  chain  rule  for  partial 
derivatives  to  achieve  this.  We  find  that 
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Substituting  for  the  quantities  g— -  ,  gx  ,  gt  9z. 

"  i 
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from  Eqs.  (3-12),  (3-13),  and  (3-19),  respectively,  we  obtain 
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where  vz  is  the  balloon  ascent  rate,  and  Wq,  W^,  and  are  defined  by  Eqs. 

(3-9),  (3-10),  and  (3-20),  respectively. 


Analogous  expressions  for  the  North  component  of  ballistic  wind  can  be 
obtained  in  the  same  fashion.  For  example, 


3V_ 
3z . 


dyi 


(5-14) 


where  Wv^  is  defined  by  Eq .  (3-24). 

When  the  values  of  the  partial  derivatives  discussed  above  are  appropriately 
substituted  into  Eq.  (5-2)  or  Eq .  (5-3),  the  resulting  expressions  are  rather 
unwieldy.  In  order  to  formulate  the  results,  we  adopt  a  mnemonic  code  to 
represent  each  of  the  18  error  suras  in  Eqs.  (5-2)  and  (5-3).  (The  first  term 
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in  brackets  in  each  equation  is  taken  to  be  a  single  error  sum.)  This  code  is 
displayed  in  Table  I.  Code  names  beginning  with  the  letter  B  represent  bias 
sums  before  squaring  and  thus  have  units  of  velocity;  code  names  beginning 
with  the  letter  R  represent  sums  of  the  squares  of  random  error  contributions 
and  thus  have  units  of  velocity  squared. 


The  error  sums,  after  substitution  for  the  partial  derivatives  of  U  and  V,  are 
displayed  in  Tables  II  and  III.  The  forms  of  the  partial  derivatives  of  the 

3x^  3x^ 

geometric  variables,  e.g.,  y—  ,  —  etc..,  are  given  in  Appendix  B. 

In  RAWIN  (and  also  in  the  RADAR  model),  it  is  assumed  that  single  values  are 
valid  over  all  zones  for  the  following:  bias  error  og^  in  elevation  due  to 
causes  other  than  ground  reflection;  random  error  Orel  *-n  elevation  due  to 
causes  other  than  ground  reflection;  bias  error  age  in  azimuth;  and  random 
error  oga  in  azimuth.  The  remaining  errors  may  vary  from  zone  to  zone  and 
hence  retain  the  subscript  i  in  Tables  II  and  III. 


Each  of  the  errors  oBelG  and  oR£ 1G  due  to  ground  reflection  depends  on  (among 
other  things)  the  angle  of  incidence  that  the  incoming  signal  makes  with 
respect  to  the  surface  of  the  earth.  If  the  terrain  in  front  of  the  receiving 
set  is  sloping,  the  angle  of  incidence  will  be  affected.  This  should  be  taken 
into  account  in  calculating  the  sums  BERXE,  RERXE,  BERYE,  and  RERYE.  In  RAWIN 
(and  also  in  RADAR),  a  positive  or  zero  angle  of  slope,  called  the  foreground 
elevation  F,  is  required  as  input.  The  model  utilizes  F  to  determine  the 
correct  angle  of  incidence.  It  then  uses  linear  interpolation  to  select 
proper  values  of  ajjeiG  ar*d  °ReiG  from  inputted  tables.  These  tables  must 
contain  externally  generated  values  of  ogc jG  anc*  aRejG»  j  =  1.  2,  •••,  271, 
corresponding  to  potential  elevation  angles  of  0.0,  0.33,  0.67,  1.0,  1.33, 
...,  90.0  degrees,  respectively. 


Table  IV  summarizes  all  the  inputs  required  by  the  RAWIN  model  for  the  case  of 
a  complete  zone  structure  containing  Nz  zones.  Conversion  of  input  units  to 
MKS  values  is  done  by  the  model  itself  where  necessary. 
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The  variance  in  the  East  component  of  ballistic  wind  is  given  Dy 


Oy2  =  REXL  +  (BEXZ)2  +  (BELXE)2  +  (BERXE)2  +  (BEXA)2 
+  REXZ  +  RELXE  +  RERXE  +  REXA. 


(5-15) 


while  the  variance  in  the  North  component  is  given  by 

Ov2  =  REYL  +  (BEYZ)2  +  (BELYE)2  +  (BERYE)2  +  (BEYA)2 

+  REYZ  +  RELYE  +  RERYE  +  REYA. 


(5-16) 


The  CW  is  computed  from  Eq.  (4-13).  On  output  the  units  are  converted  to 
knots  squared. 


Table  I.  Mnemonic  Code  for  RAWIN  Error  Sums 


East  North 


Error  Source 

Bias 

(m/ sec) 

Random 

(m/sec)2 

Bi  as 

( m/sec ) 

Random 

(m/sec)' 

Launch  Position  Determination 

-- 

REXL 

-- 

REYL 

Altitude  Determination 

BEXZ 

REXZ 

BEYZ 

REYZ 

Elevation  Measurement  ("Leveling") 

BELXE 

RELXE 

BELYE 

RELYE 

Elevation  Measurement  (Reflection) 

BERXE 

RERXE 

BERYE 

RERYE 

Azimuth  Measurement 

BEXA 

REXA 

BEYA 

REYA 
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Table  II. 


RAWIN  Error  Sums  for  the  East  Component  of  Ballistic  Wind  (Line  n) 
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Table  III.  RAWIN  Error  Sums  for  the  North  Component  of  Ballistic  Wind  (Line  n) 
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1  z  0  3D„ 


2  3yn  2 

*  °RD0  +  ^  V~W'>  '  a° 


z  0  daQ  '  RaO 


n  3y . 

BEYZ  =  Z  [  v  W.  ^  +  W  .  ]  an  . 

z  i  3z .  vi  Bzi 


i=l 


REYZ  =  Z  [  v  W. 
i=l 


3y; 


+  w.._.  I  a* 


z  i  3z .  vi  1  Rzi 

i 


n  3y . 

BELYE  =  Z  v  W.  5-^  an  T 

.  ,  z  i  3e .  BeL 

i  =  l  i 


n  3yi  2  o 

REI.YE  =  2  (  vzW.  1 

1  =  1  l 


n  3y 

BERYE  =  Z  v  W.  g-i  o„  .  r 
z  i  3e .  BeiG 
i=l  i 


RERYE  = 


n  3y .  2 

1  [  vzWi  3e7  *  CTReiG 

i  =  l  i 


n  ■  3y . 

BEYA  =  I  v  W.  ^  a, 


i=l 


z  i  3a  Ba 
1 


n  3y  2 

REyA  =  Z  [  v  W  5-  ]  cr  2 


i=l 


z  i  da .  1  Ra 

i 
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Table  IV.  RAWIN  Inputs 


Symbol 

Input 

Units 

No .  of 
Values 

Explanation 

Do 

m 

1 

Launch  Displ.:  used  to  compute  partials 

F 

deg 

1 

Foreground  Elev. :  used  in  selecting  values  of 
°BeiG  and  °ReiG  from  aBejG  aad  aRejG 

u . 

1 

m/sec 

N 

z 

E  st  zone  wind  vel.:  used  to  compute  W  ^ 

V  . 

1 

m/sec 

N 

z 

North  zone  wind  vel.:  used  to  compute  W 

V 

z 

m/min 

1  to  4 

Balloon  ascent  rate 

wik 

N  2 
z 

Complete  table  of  wind  weighting  factors,  where 
i  =  1,  ...,  and  k  =  I,  ...,  ;  for  any  line  n 

the  model  selects  approp.  values  of  w.  which  are 

used  to  compute  W. ,  W  ,  W  . ,  W  . . 

r  O’  i’  ui  vx 

z . 

1 

m 

N 

z 

Zone  top  alt.:  used  to  compute  W. ,  W.,  W  W  . 

r  r  0’  i’  ui  vi 

ao 

deg 

1 

Launch  azimuth:  used  to  compute  partials 

a . 

i 

deg 

N 

z 

Zone  top  azimuth:  Used  to  compute  partials 

£  . 

1 

deg 

N 

z 

Zone  top  elev.:  used  to  compute  partials 

ctrdo 

m 

1 

Random  error  in  launch  displacement 

°RaO 

deg 

1 

Random  error  in  launch  azimuth 

°Bzi 

m 

N 

z 

Bias  error  in  altitude 

CTRzi 

m 

N 

z 

Random  error  in  altitude 

aBa 

deg 

1 

Bias  error  in  azim.  tracking  of  apparent  target 

aRa 

deg 

1 

Random  error  in  azim.  tracking  of  apparent  target 

°BcL 

deg 

1 

Bias  error  in  elev.  tracking  of  apparent  target 

aReL 

deg 

1 

Random  error  in  elev.  tracking  of  apparent  target 

°BejG 

deg 
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Potential  bias  errors  in  elev.  due  to  ground 
reflection:  used  to  determine 

°RejG 

deg 

271 

Potential  Random  errors  in  elev.  due  to  ground 
reflection:  used  to  determine  a_ 

K£  10 

NOTE :  N 

=  15  for 

the  NATO 

zone  structure. 

6.  RADAR 


RADAR  models  the  case  in  which  a  ground-based  radar  set  measures  the  indepen¬ 
dent  variables  of  slant  range,  elevation  angle,  and  azimuth  associated  with 
the  radiosonde’s  position.  This  model  is  applicable,  for  example,  to  the 
Radar  Sets  AN/TMQ-19  and  AN/FPS-16. 

The  derivation  of  the  variance  in  each  component  of  ballistic  wind  follows 
essentially  the  same  path  as  in  RAWIN.  The  main  difference  in  RADAR  is  that 
the  altitude  z^  at  zone  top  i  is  not  considered  to  be  an  independent  variable 
but  rather  is  a  function  of  slant  range  St  and  the  elevation  angle  e^.  Then 
from  Eq.  (5-1),  we  may  indicate  the  functional  form  of  the  East  component  U  of 
ballistic  wind  by 


U  =  U  [xQ  (D, 


a0). 


c.  (z.  (S.,  6.),  e.,  a.),  t.  (z.  (S., 


Ei)) 


(6-1) 


i=  1 >  2 ,  ...»  n  , 


where  all  symbols  retain  their  previously  defined  meanings. 

It  is,  of  course,  not  necessary  to  write  U  in  precisely  the  form  shown  in  Eq . 
(6-1).  We  could,  for  example,  omit  the  intermediate  variable  zj  entltely. 
However,  Eq.  (6-1)  permits  us  to  make  further  use  in  radar  of  several  partial 
derivatives  that  are  also  valid  in  RAWIN. 

The  variance  o^2  in  U  is  obtained  by  specializing  Eq.  (4-3)  to  the  independent 
variables  D0,  oQ,  S±,  and  cy  The  result  is 


■(§U_0  )2 

■l3D0  RD0J 


,au  s2^ 
(daQ  aRa (P  1 


+ 
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n 

I  I 

i=l 


au_ 

dS . 


BSi 


n 

♦  1  * 
i=l 


au 

dc 

i 


BeiL1 


n 

(  2 
i  =  l 


au 

dc 

1 


BeiG 


I2  +  i 


n 

2 

i  =  l 


au 

3a 


Bai 


(6-2) 


n 

2 

i=l 


[3U  „  , 

Us.  CTRSi 
1 


a 

2 

i=l 


,3U  „  12  + 

3e .  °ReiL 
1 


n 

I 

i=l 


&  °R£iGl2  + 

1 


a 

I 

i=  1 


[§U-  0  ] 2 

l3a.  Rai1 
i 


In  Eq.  (6-2),  each  type  of  error  is  designated  by  a  with  appropriate  sub¬ 
scripts.  The  meanings  of  the  subscripts  are  identified  by  the  rules  that  are 
stated  below  Eq .  (5-2).  For  example,  oRS^  is  the  random  error  in  the  measure¬ 
ment  of  slant  range  for  the  i'th  zone  top.  As  before,  we  assume  that  any 
error  in  the  launch  time  is  negligible  and  omit  a  term  involving  to  in  Eq. 
(6-2). 

Similarly,  the  variance  av2  in  th.'  North  component  V  of  ballistic  wind  is 


CTV2  =  t(3D0  °RD0)2  +  (aa0  °Ra0^2 ^ 


+  1  I 


av 


,  as.  BSi 

1=1  i 


a_.l2  +  [  2 


av 


.  ,  dc.  BeiL 
i=l  i 


2  +  [  2 


i=l  i 


liT  W*  +  (  * 


i=l  i 


o  ]2 

3a.  BaiJ 


(6-3) 


n  ,8V  _  ,,  .  ?  r9V 


*  *  i&rW  VsrW* 

i=l  i  i=l  i  i-i  i  1-1 


,3V  )2 

3a .  aRai 
1 


We  use  the  chain  rule  for  partial  derivatives  in  conjunction  with  Eq.  (6-1)  to 
obtain 


au 

3x 


^0 
o  aDo 


(6-4) 
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(6-5; 


9U 

9x, 


3x, 


3a, 


3U 

3U 

3x . 

1 

3z . 

1 

3U 

3t . 
_ 1 

3z . 

1 

3S . 

1 

3x . 

1 

3z . 

1 

3S . 

1 

3z . 

1 

3S . 

1 

> 

3U 

3U 

3x . 

1 

3z. 

1 

+  3U_ 

3x . 

1 

+  §3L 

3t .  3z . 
_ 1  _ 1 

3e . 

1 

3x . 

L 

3z . 

1 

3c. 

1 

3x . 

1 

3c. 

1 

3t . 

1 

3z .  3e . 

1  1 

and 


3U 

3a 


9U 
9x . 


With  appropriate  substitution  from  Eqs .  (3-12),  (3-13),  and  (3-19),  the 
expressions  become 


3U_  % 

3D0  =  VZW0  3Dq 


3U  _  M  3X0 
3a0  -  Vz  0  3a0 


3U 
3S . 


3x .  3z . 

(v  W.  5—^  +  W  . ) 

z  1  3z.  ui  3S . 

1  1 


3U 

3e . 
1 


3x .  3z .  3x . 

(v  W.  —  +  W  .)  +  v  W. 

z  1  3z .  ui  3e.  -  ’ 

1  1 


z  i  3e . 

1 


and 


3U 

3a. 


3x . 

v  W.  x — - 
z  1  3a . 

1 
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(6-6) 


(6-7) 


(6-8) 


above 


(6-9) 


(6-10) 


(6-11) 


(6-12) 


(6-13) 


where  vz  is  the  balloon  ascent  rate,  and  Wo,  ,  and  Wuj  are  defined  by  Eqs. 
(3-9),  (3-10),  and  (3-20),  respectively. 

Analogous  expressions  involving  the  North  component  of  ballistic  wind  can  be 
obtained.  For  example. 


9V 

3S 


3yi 


(6-14) 


where  Wv^  is  defined  by  Eq .  (3-24). 

The  18  error  sums  in  Eqs.  (6-2)  and  (6-3)  are  represented  by  the  mnemonic  code 
presented  in  Table  V.  The  algebraic  forms  of  these  sums,  after  substitution 
for  the  various  partial  derivatives  of  U  and  V,  are  given  in  Tables  VI  and 
VII. 

Since  the  altitudes  2|  define  the  zone  structure,  they  are  still  required  as 
inputs  to  RADAR.  Conversely,  actual  values  of  slant  range  are  not 
required.  Such  values  could  be  used  to  compute  some  of  the  partial  deriva¬ 
tives  shown  in  Tables  VI  and  VII.  However,  values  of  zj  also  serve  in  this 
capacity  and  are  used  instead.  See  Appendix  B. 

The  discussion  near  the  end  of  Section  5  concerning  errors  in  azimuth,  errors 
in  elevation,  and  the  foreground  elevation  is  also  pertinent  to  the  RADAR 
model.  In  addition,  RADAR  requires  a  single  inputted  value,  valid  for  all 
zones,  for  the  random  error  Ogg  in  slant  range  and  also  a  single  value  for  the 
bias  error  ORS*  Although  it  might  be  expected  that  ogg  would  be  negligible, 
it  is  nevertheless  included  here  for  generality. 

Table  VIII  summarizes  all  the  inputs  required  by  the  RADAR  model  for  the  case 
of  a  complete  zone  structure  containing  Nz  zones.  Conversion  of  input  units 
to  MRS  values  is  done  by  the  model  itself  where  necessary. 
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The  variance  in  the  East  component  of  ballistic  wind  is  given  by 


ou2  =  RE XL  +  ( BEXS 1 2  +  (BELXE)2  +  (BERXE)2  +  (BEXA)2 

+  REXS  +  RELXE  +  RERXE  +  REXA  , 

while  the  variance  in  the  North  component  is 

Oy2  =  REYL  +  (BEYS)2  +  (BELYE)2  +  (BERYE)2  +  (BEYA)2 

+  REYS  +  RELYE  +  RERYE  +  REYA 

The  CW  is  computed  from  Eq.  (4-13).  Output  is  in  knots  squared. 


(6-15) 


(6-16) 


Table  V.  Mnemonic  Code  for  RADAR  Error  Sums 


East  North 


Error  Source 

Bias 

(m/sec) 

Random 

(m/sec)2 

Bias 

(m/sec) 

Random 

(m/sec): 

Launch  Position  Determination 

— 

REXL 

-- 

REYL 

Slant  Range  Measurement 

BEXS 

REXS 

BEYS 

REYS 

Elevation  Measurement  ("Leveling") 

BELXE 

RELXE 

BELYE 

RELYE 

Elevation  Measurement  (Reflection) 

BERXE 

RERXE 

BERYE 

RERYE 

Azimuth  Measurement 

BEXA 

REXA 

BEYA 

REYA 
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Table  VI.  RADAR  Error  Sums  for  the  East  Component  of  Ballistic  Wind  (Line  n) 


RE  XL 

= 

[  v 

3xo 

zW°  3Dq  ' 

2 

1  a 

2 

RDO 

+  [ 

O 

N 

> 

8xo 

2 

] 

a  ^ 

RaO 

n 

9x . 

9z . 

BEXS 

2 

i=l 

[( 

1 

82 . 

1 

+ 

W  .  ) 

Ul 

1 

3S  . 

1 

)  < 

3bs 

n 

9x . 

9z  . 

2 

REXS 

2 

i  =  l 

[(  v^. 

1 

82 . 

1 

+ 

W  .  ) 

Ul 

1 

9S . 

1 

] 

°RS 

n 

9x . 

9z  . 

3x . 

BELXE 

2 

i=l 

[  (  vzW. 

1 

32. 

1 

+ 

W  .  ) 

Ul 

1 

9e . 

1 

+ 

v  W . 
z  1 

1 

9c  . 

1 

1  °BeL 

n 

3x  . 

9z . 

9x . 

2 

RELXE 

2 

i=l 

I(  VzWi 

1 

82 

1 

+ 

W  .  ) 

Ul 

1 

3c. 

1 

+ 

v  W. 
z  1 

1 

3c . 

1 

1  °RcL 

n 

3x . 

9z . 

9x . 

BERXE 

2 
i  =  1 

[( 

1 

92. 

1 

w  .  ) 

Ul 

1 

3c. 

1 

+ 

v  W . 
z  1 

1 

9e" 

1 

1  °B£iG 

n 

3x . 

9z . 

Bx . 

2 

RERXE 

2 

i=  1 

((  v^. 

1 

9e . 

1 

+ 

WU1  } 

1 

9 FT 

1 

+ 

v  W . 
z  1 

1 

3c. 

1 

1  °ReiG 

n 

3x 

BEXA 

I 

i=l 

v  W .  5- 

z  1  3a 

i 

-  ( 

i 

JBa 

n 

3x . 

z 

REXA 

- 

2 
i=  1 

[  v  W . 
z  1 

1 

9a . 

1 

1 

°Ra 
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Table  VII. 


RADAR  Error  Sums  for  the  North  Component  of  Ballistic  Wind  (Line  n) 


REYL  = 

1  vz 

3y0  2 

W  — -  ) 

0  3Dq 

a  2  + 

RDO 

1  VzW0 

9yo 

9ao 

'  °RaO 

n 

3yi 

3z  . 

i  , 

2 

BEYS  = 

2 

i  =  l 

[  (  v  W. 

2  1 

— -  + 
3z . 

l 

Wvi  > 

3S.  J 
i 

°BS 

REYS  = 

n 

2 

i  =  l 

[ (  VzWL 

^  ♦ 
3z 

i 

w  .  ) 

VI 

3z . 

—L  ) 

3S.  J 
i 

°Ri 

n 

3yi 

3z  . 

i 

+  v  W. 

Z  1 

9yL 

BELYE  = 

2 

i=l 

(( 

— i  + 

3z  . 

l 

w  .  ) 

VI 

3&T 

l 

3e . 

i 

n 

8yi 

3z . 

+  v  W. 

Z  1 

9yi 

RELYE  = 

I 

i=l 

[(  vzW. 

3z7  + 

i 

W  .  ) 

VI 

3 T 

3e . 

i 

n  aVi  azi  3yi  , 

BERYE  =  2  [(  vzW.  gj-  +  Wv.  )  +  vzWi  9e .  1  ViG 

i=l  i  11 


n 

3yi 

3z . 

i 

u  9yi 

RERYE  = 

2 

i=l 

“  VzWi  8eT  + 

1 

w  .  ) 

VI 

3FT  + 

i 

vzwi  §r 

i 

n 

9yi 

BEYA  = 

2 

i=l 

vzWi  35“  aBa 
1 

REYA  = 

n 

2 

i=l 

2 

1  Vi  SJ-> 

1 

aRa 
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Table  VIII.  RADAR  Inputs 


Input 

No .  of 

Symbol 

Units 

Values 

Explanation 

Do 

m 

1 

Launch  Displ.:  used  to  compute  partials 

F 

deg 

1 

Foreground  Elev.:  used  in  selecting  values  of 

0o  . „  and  on  .  „  from  o.  ..  and  a_ 

BeiG  RtiG  B£jG  RejG 

u. 

1 

m/sec 

N 

2 

East  2one  wind  vel.:  used  to  compute 

V  . 

1 

m/sec 

N 

2 

North  2one  wind  vel.:  used  to  compute  W 

V 

z 

m/min 

1  C.0  ^ 

Balloon  ascent  rate 

w. , 
xk 

N  ‘ 

2 

Complete  table  of  wind  weighting  factors,  where 

i  =  1,  and  k  =  1,  ;  for  any  line  n 

the  model  selects  approp.  values  of  w.  which  are 

used  to  compute  W„ ,  W. ,  W  . ,  W  . .  1 

0  l  ui  vi 

z . 

X 

m 

N 

2 

Zone  top  alt.:  used  to  compute  to’  ,  to.  ,  to  .  ,  to  . 

ao 

deg 

1 

Launch  azimuth:  used  to  compute  partials 

a . 

i 

deg 

N 

2 

Zone  top  azimuth:  used  to  compute  partials 

£. 

1 

deg 

N 

2 

Zone  top  elev. :  used  to  compute  partials 

°RD0 

m 

1 

Random  error  in  launch  displacement 

aRaO 

deg 

1 

Random  error  in  launch  azimuth 

°BS 

m 

1 

Bias  error  in  slant  range 

aRS 

m 

1 

Random  error  in  slant  range 

aBa 

deg 

1 

Bias  error  in  azim.  tracking  of  apparent  target 

JR  a 

deg 

1 

Random  error  in  azim.  tracking  of  apparent  target 

aB£L 

deg 

1 

Bias  error  in  elev.  tracking  of  apparent  target 

CTR£L 

deg 

1 

Random  error  in  elev.  tracking  of  apparent  target 

aB£jG 

deg 

271 

Potential  bias  errors  in  elev.  due  to  ground 

reflection:  used  to  determine  a„ 

BeiG 

aRejG 

deg 

271 

Potential  Random  errors  in  elev.  due  to  ground 
reflection:  used  to  determine  CT_  . _ 

K£lu 

MOTE :  N 

z 

=  15  for 

the  MATO 

zone  structure. 
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7 .  NAVAID 


NAVAID  models  the  case  In  which  radionavigation  techniques  are  used  to  deter¬ 
mine  the  position  of  the  radiosonde.  See  Reference  3,  for  example,  for  a 
review  of  these  techniques.  NAVAID  is  a  simplified  model  in  that  it  bypasses 
the  complexities  of  hyperbolic  geometry  by  requiring  estimated  errors  in  the 
East  and  North  components  of  position  as  input.  This  model  is  applicable  to 
the  radionavigation  portions  (LORAN,  VLF,  and  OMEGA)  of  FAMAS. 

To  obtain  the  implicit  functional  form  of  the  East  component  U  of  ballistic 
wind  appropriate  to  line  n,  we  rewrite  Eq.  (4— 1 )  in  the  following  way: 


U  U  (xq>  t0’  xi ’  (z^)]>  i  =  2, 


(7-1) 


The  independent  variables  are  taken  to  be  the  East  launch  coordinate  xq,  the 
launch  time  to,  the  radiosonde's  East  coordinate  x^  at  the  top  of  zone  i,  and 
the  corresponding  altitude  z^.  As  in  RAWIN,  the  time  t^  is  considered  to  be  a 
function  of  z^. 

In  hyperbolic  tracking  any  bias  errors  in  and  in  y^,  the  North  component  of 
position,  are  expected  to  be  negligible.  Therefore,  bias  error  sums  are  not 
computed  for  these  variables  in  NAVAID.  However,  we  do  allow  for  both  bias 
and  random  errors  in  the  determination  of  z^.  As  in  the  other  models,  it  is 
assumed  that  there  is  no  error  associated  with  tg* 

Specializing  Eq.  (4-3)  to  NAVAID,  we  have  for  the  variance  o^2  in  U: 


2  - 


=  [ 


3U 

3x 


0 


RxO 


r  v  9U 

*  *  Bz.  °Bzi 

i=l  i 


n 

+  I 
i  =  l 


n 

I 

i=l 


BU 
Bz . 


Rzi 


(7-2) 
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where  each  type  of  error  is  designated  by  a  with  appropriate  subscripts.  The 
meanings  of  the  subscripts  are  identified  by  the  rules  stated  below  Eq. 
(5-2).  For  example,  Is  the  random  error  in  the  East  coordinate  for  zone 
top  i.  A  similar  expression  holds  for  the  variance  av2  in  the  North  component 
V  of  ballistic  wind: 


3 V_ 
9y0 


°RyO 


+ 


n 

I 

i=l 


n 

+  I 
i=  1 


av 

9y  Ryi 


n 

+  I 
i=l 


av 

3z . 


Rzi 


(7-3) 


The  mnemonic  code  representing  the  eight  error  sums  in  Eqs.  (7-2)  and  (7-3)  is 
given  in  Table  IX. 


The  partial  derivatives 


are  given  by  Eqs.  (3-12), 


au  au  av  av 

ax0  *  a*,  •  ay0  ’  an  ay, 

(3-13),  (3-21),  and  (3-22),  respectively.  Using  the  chain  rule  for  partial 
derivatives  and  substituting  from  Eqs.  (3-19)  and  (3-23),  we  also  have 


au 

3z. 


w  . 

Ul 


(7-4) 


and 


av 

3z . 


W  . 
vi 


(7-5) 


where  and  Wv,  are  defined  by  Eqs.  (3-20)  and  (3-24),  respectively. 

NAVAID  computes  the  random  positional  errors  ciRxi  and  using  th  .ocedure 

outlined  below. 
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It  is  assumed  that  the  hyperbolic  tracking  system  fixes,  i.e.,  determines,  the 
East  and  North  components  of  the  radiosonde's  position  every  Ax  seconds.  We 
take  the  random  error  in  the  fixing  of  each  component  to  be  a  constant  that  is 
characteristic  of  the  measuring  system.  These  errors,  labeled  <jrxt  and  ORy-j , 
respectively,  are  required  as  inputs  to  NAVAID. 

Consider  a  height  Interval  of  AH^  meters,  which  is  centered  at  the  top  of  zone 
i.  As  the  radiosonde  traverses  AH^,  fixes  of  position  are  made.  We  have 


N 


i 


AH. 

Integer  [  ^  ] 
z 


(7-6) 


where  the  ascent  rate  vz  is  expressed  in  meters  per  second. 

The  value  of  x^  at  the  zone  top  is  taken  to  be  the  mean  of  the  N^  fixes  of  the 
radiosonde's  East  coordinate.  Under  this  condition  is  given  by 


Rxi 


Rxx 

N. 


Similarly, 


Ry  i 


Ryt 


(7-7) 


(7-8) 


The  component  launch  errors  o^q  anc*  aRyO  are  I130^6^  *n  tI,e  following  way. 
As  the  radiosonde  sits  on  the  ground,  its  position  can  be  determined  by  hyper¬ 
bolic  fixing  during  some  time  interval.  The  radiosonde-balloon  combination 
may  then  be  transported  to  another  nearby  position  for  actual  launch.  The 
second  position  is  determined  relative  to  the  first  by  some  direct  method. 
Since  the  errors  in  these  two  different  types  of  measurements  are  uncorre¬ 
lated,  we  have  for  the  East  component: 
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(7-9) 


RxO 


=  or, 


RxOF 


RxOD 


where  CrxOF  is  the  component  random  error  in  the  hyperbolic  fixing  of  the 
first  position,  and  orxQD  *s  the  component  random  error  in  the  direct  deter¬ 
mination  of  the  second  position  relative  to  the  first. 

If  the  fixing  of  the  initial  position  occurs  over  a  time  interval  of  AT  min¬ 
utes,  then  the  number  Nq  of  fixes  is  given  by 

NQ  =  Integer  (60  ~]  ,  (7-10) 


and  the  variance  is 


Rxt 


RxOF 


(7-11) 


We  have  no  foreknowledge  of  the  value  of  Orx0D  nor  the  corresponding  North 
error  CgyOD ’  For  any  particular  launch  they  are  not  necessarily  equal  to  each 
other.  However,  on  average  over  many  launches,  we  expect  the  following  rela- 
t ionship: 


O 

CTRxOD 


RyOD 


(7-12) 


where  Ordq  is  the  random  error  in  the  direct  measurement  of  the  distance  from 
the  initial  position  to  the  final  launch  position.  For  the  purpose  of  model¬ 
ing  the  data  acquisition  system,  Eq .  (7-12)  is  assumed  to  hold.  The  quantity 
ordo  is  referred  to  as  the  random  launch  error  and  is  a  required  input  to 
NAVAID. 
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From  Eqs .  (7-9),  (7-11),  and  (7-12),  we  have 


RxO 


o  2  a  2 

Rxt  RDO 

2 


An  analogous  equation  yields  o^y 


RyO. 


(7-13) 


In  order  to  determine  N^  and  Nq  ,  NAVAID  takes  At  =  1  second  and  AT  =  5  min¬ 
utes.  It  also  assumes  the  following  arbitrary  values:  AHj  =  200  meters  for 
the  first  five  zones;  AHj  =  400  meters  for  the  remaining  zones.  The  motiva¬ 
tion  for  this  is  the  NATO  zone  structure  in  which  the  higher  zones  are  signif¬ 
icantly  thicker  than  the  lower  zones.  The  Users'  Manual  explains  how  to  alter 
the  values  of  At,  AT,  and  AH^. 

The  preceding  analysis  is  used  to  find  the  final  form  of  each  of  the  required 
bias  and  random  sums.  These  are  displayed  in  Table  X. 

The  inputs  to  NAVAID  are  summarized  in  Table  XI.  Conversion  of  input  units  to 
MKS  values  is  done  by  the  model  itself  where  necessary. 

The  variances  in  the  East  and  North  components,  respectively,  of  ballistic 
wind  are  given  by 


dy2  =  REXL  +  (BEXZ) 2  +  RE XX  +  REXZ 


(7-14) 


and 


ov2  =  REYL  +  (BEYZ)2  +  REYY  +  REYZ 


(7-15) 


The  CW  is  computed  from  Eq.  (4-13).  Output  is  in  knots  squared. 


Table  IX.  Mnemonic  Code  for  NAVAID  Error  Sums 


East 


North 


Error  Source 


Bias  Random  Bias  Random 

(m/sec)  (m/sec)2  (m/sec)  (m/sec)2 


Launch  Position  Determination 


REXL  --  REYL 


Altitude  Determination 


BEXZ 


REXZ 


BEYZ 


REYZ 


r 


Table  XI.  NAVAID  Inputs 


t 

i 


Symbol 

Input 

Units 

No .  of 
Values 

Explanation 

u . 

1 

m/sec 

N 

z 

East  zone  wind  vel.:  used  to  compute 

V  . 

1 

ra/sec 

N 

z 

North  zone  wind  vel.:  used  to  compute  W 

V 

z 

ra/min 

1  to  4 

Balloon  ascent  rate 

w. . 
ik 

N 

z 

Complete  table  of  wind  weighting  factors,  where 

i  =  1,  .  . . ,  N^  and  k  =  1,  ....  N^;  for  any  line 

the  model  selects  approp.  values  of  w^ ,  which 

are  used  to  compute  W„ ,  W. ,  W  W  , . 

r  O’  i*  ui’  vi 

z . 

t 

m 

N 

z 

Zone  top  alt.:  used  to  compute  W„,  W,  W  W  . 

r  0’  i’  ui ’  vi 

aRDO 

m 

1 

Random  error  in  direct  determination  of  launch 
position 

QBzi 

m 

N 

z 

Bias  error  in  altitude 

°Rzi 

m 

N 

z 

Random  error  in  altitude 

aRxt 

m 

1 

Random  error  associated  with  a  single  hyperbolic 
fix  of  the  East  coordinate 

aRyt 

m 

1 

Random  error  associated  with  a  single  hyperbolic 
fix  of  the  North  coordinate 

NOTE :  N 

=  15  for 

the  NATO 

zone  structure. 

z 
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8.  COMMENTS  AND  SUGGESTIONS  FOR  FURTHER  WORK 


In  RAWIN  and  RADAR,  errors  in  elevation  due  to  ground  reflection  can  make 
significant  contributions  to  the  error  in  ballistic  wind,  particularly  for  the 
case  of  low  elevation  angles.  These  errors  depend  on  a  number  of  factors, 
including  elevation  angle,  surface  dielectric  constant,  signal  frequency, 
antenna  voltage  pattern,  and  the  tracking  method  used.  General  information  on 
two  tracking  methods,  sequential  lobing  and  conical  scan,  can  be  found  in 
Reference  4.  A  limited  auxiliary  model  LRDC  exists  for  computing  theoretical 
values  of  bias  and  random  reflection  errors  for  an  antenna  pattern  that  is 
indentical  for  the  high  and  low  switched  positions;  its  implementation  is 
described  in  the  Users'  Manual. 

The  bias  reflection  error  generated  by  LRDC  for  any  given  elevation  angle  does 
not  represent  a  standard  deviation  characteristic  of  a  normal  population  with 
zero  mean  error.  Consequently,  if  LRDC  is  used  to  generate  the  required  bias 
reflection  inputs  to  RAWIN  and  RADAR,  it  should  be  recognized  that  the  quanti¬ 
ties  o^2  and  Oy 2  discussed  in  this  report  no  longer  correspond  precisely  to 
variances.  (The  bias  errors  computed  by  LRDC  are  typically  much  smaller  than 
the  corresponding  random  reflection  errors,  which  are  also  computed.) 

Most  of  the  error  sums  obtained  in  this  report  for  ballistic  line  n  contain 
the  factor  W^,  1  =  1,  2,  ...,  n.  As  defined  by  Eq .  (3-10),  the  quantity 
represents  the  wind  weighting  factor  per  unit  zone  width  for  zone  i  minus  the 
wind  weighting  factor  per  unit  zone  width  for  zone  l  +  1.  Depending  on  the 
weighting  factors  and  zone  widths  involved,  W^  can  be  positive,  negative,  or 
zero.  For  the  case  of  Message  3  weighting  factors,  Wn  has  the  largest  magni¬ 
tude  of  any  member  of  the  set  {W^,  i  =  1,  2 .  n}  for  the  following 

reasons:  the  zone  wind  weighting  factor  is  greatest  for  zone  n;  the  wind 
weighting  factor  for  zone  n  +  1  is  zero.  For  example,  using  the  NATO  zone 
structure  and  the  corresponding  Message  3  wind  weighting  factors  given  in 
Reference  2,  one  obtains  the  following  values  (in  meter  *)  relevant  to  ballis¬ 
tic  line  7:  »l  -  0.;  W2  -  -0.40  x  10~4;  W3  =*  0.;  W 4  -  -0.20  x  10-4;  W&  = 
-0.40  x  10  4;  Wc  =  -0.33  x  10  W,  ■  0.53  x  10  The  way  in  which  the  W. 
are  calculated  results  in  a  "cancellation  effect"  for  i  less  than  n.  In  fact, 
if  Wj  *  0.,  the  contribution  of  zone  i  to  some  of  the  error  sums  will  be 
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nullified.  This  is  not  necessarily  what  one  would  expect  intuitively.  In 
order  to  verify  whether  or  not  the  cancellation  effect  is  real,  appropriate 
comparisons  between  model  predictions  and  experimental  results  should  be 
attempted. 

Each  of  the  models  utilizes  a  single  zone  wind  profile,  u^,  v^,  1  *  1,  2 . 

N2  (where  nz  is  the  total  number  of  zones  in  the  structure),  valid  for  all 
balloon  ascent  rates.  From  a  known  or  postulated  profile,  tables  of  values  of 
elevation  and  a^,  i  *  1,  2,  ...,  Nz  can  be  developed  fo.  each  balloon 
ascent  rate  for  use  in  RAW1N  or  RADAR.  Alternatively,  if  the  values  of  Ej  and 
are  known,  a  zone  wind  profile  can  be  calculated.  In  either  case,  the 
tables  of  values  must  be  obtained  externally  to  the  model  and  used  as  input  to 
it.  Although  the  procedure  for  doing  this  is  conceptually  straightforward,  it 
would  be  useful  to  incorporate  it  into  the  models  themselves  or,  alterna¬ 
tively,  create  an  auxiliary  model  that  would  generate  the  necessary  tables. 

The  zone  wind  profile  is,  of  course,  merely  a  useful  representation  of  actual 
wind  conditions  aloft.  Real  winds  do  not  necessarily  maintain  constant  magni¬ 
tude  and  direction  over  an  arbitrary  zone  width  and  then  abruptly  assume  new 
values  at  a  zone  boundary.  Thus  the  zone  wind  methodology  Introduces  an 
artificial  discontinuity  or  "shear"  in  wind  at  the  zone  boundaries.  Any  error 
in  the  determination  of  the  zone  top  altitudes  will  lead  to  a  calculated  zone 
wind  profile  that  does  not  correspond  precisely  to  the  zone  structure  under 
study.  The  form  of  the  quantities  Wuj  and  Wvj  in  Eqs.  (3-20)  and  (3-24), 
respectively,  and  the  manner  in  which  they  were  obtained  suggest  that  they  are 
related  to  first-order  wind  shear  contributions  to  the  error  in  ballistic 
wind.  Further  work  on  this  point  and  on  the  entire  question  of  wind  shear  is 
suggested. 

The  models  presented  here  represent  one  approach  to  the  subject  of  ballistic 
wind  measurement  error  analysis.  They  are  not,  of  course,  the  final  w>rd  on 
this  subject.  The  subsequent  use  of  these  models  should  help  to  establish 
their  strong  points,  as  well  as  locate  areas  in  need  of  improvement  or  exten¬ 
sion. 
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APPENDIX  A.  VARIANCE  ASSOCIATED  WITH  A  COMPONENT  OF  BALLISTIC  WIND 


Using  appropriate  assumptions,  we  want  to  show  that  Eq.  (4-3)  yields  the 
variance  in  the  East  component  U  of  ballistic  wind  for  line  n. 

The  generalized  implicit  functional  form  of  U  is  given  by  Eq .  (4-2): 


U  -  U  fAlO’  X20' 


^ 1 i ’  ^2i ’ 


i-1 >  2, 


The  subsequent  analysis  is  simplified  somewhat  if  the  number  of  subscripts  in 
Eq.  (A-l)  is  reduced.  To  this  end  we  rewrite  the  subset  of  independent  launch 
variables  X^q,  ^20  ^L0>  as  V0>  Vq  ,  •••»  and  we  rewrite  the  subset  of 
independent  variables  £  ,  •••»  determined  for  the  top  of  zone  i,  as 
Pi,  Ci*  •••,  explicitly  retaining  only  two  members  of  each  subset.  Then, 


=  U  (  Mr 


=  1,  2, 


(A-2) 


In  order  to  obtain  the  expression  for  the  variance  in  U,  we  consider  a 
Gedanken  experiment  in  which  a  large  number  M  of  balloon  flights  is  carried 
out  under  the  same  meteorological  conditions.  For  the  m’th  flight  the  inde¬ 
pendent  variables  listed  in  Eq .  (A-2)  have  the  measured  values  POin,  v0m»  •••» 
pim>  Cim>  •••>  respectively.  These  can  be  used  to  calculate  a  component  value 
Um  of  ballistic  wind,  appropriate  to  flight  m. 

We  associate  U  with  the  mean  of  the  set  of  values  Um,  m=l,  2,  ...,M.  Let 
6Um  be  the  deviation  of  Um  from  the  mean.  To  first  order  <5Um  is  given  by  the 
Taylor  expansion: 
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,,,  _  au  ,  __  au  . 

6um  -  3^  6p0m  +  6vOm  + 


"  au  ,  ~  au  -y 

+  *  5—  6p.  +  I  5f~  6£.  + 

.  ,  3p.  im  ,  9t .  im 

i=l  i=l  l 


(A-3) 


where  each  of  the  quantities  6u„  ,  6v.  ,  . . . ,  6p .  ,  6t .  ,  ...  represents  the 

n  'um’  Om  Kira’  '’im 

deviation  from  the  mean  in  the  corresponding  variable.  The  partial  derivatives 
are  evaluated  using  the  mean  value  of  each  variable  appropriate  to  zone  top  i 
or  to  launch  position. 

For  large  M  the  variance  a^2  in  U  is  equivalent  to  the  mean  squared  deviation: 


=  H  I  (6U  )2 


m=l 


Thus  , 


(A-4) 


2  _ 


2  "  .  |  |!L  4,.  .  S’  i,„  * 

M  .  0pA  Om  3v_  Om 
m=l  0  0 


*  T  3U  +  v  3U  Kf  *  12 

+  I  = - Op .  +  I  ^zr-  OC  .  +  .  .  . 

9p.  im  .  .  0c  im 
i=l  'i  i=l  *i 


Expansion  of  the  square  in  Eq.  (A-5)  yields 


(A-5) 


1  M  911  1  M 

1  i  [  ^  6Mnm  ]2  +  h  l 


M  m=l  1  3p0  r°m 


M 


m=l 


jp-  6v  ]2  + 

9v0  Om 


,  M  n  .  M 

+  J  I  [I  fp-  6p.  ]2  +  £  1 

M  m=l  i=l  3p  im  M  m=l 


*  IT  «i.  + 

i  =  l  *i 


+  CROSS  TERMS 


(A-b) 
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One  of  the  cross  terms  in  Eq.  (A-6)  is  given  by  C^,  where 


r  -  2  ”  r  y  5  au. 

Cp£  M  1  f  1  3p .  <5pim  l  3  £.  6^im 

m=l  i=l  'i  i=l  bi 


Expansion  of  the  sums  in  Eq.  (A-7)  yields 


r  _  2  ”  ,  "  "  au.  au 

CpC  M  2  [  1  Z  3p .  3T  6Pim  6^jm 

m=l  i=l  j=l  Ki  J 


For  this  form  the  sum  over  m  may  be  performed  first: 


r  -  o  y  ?  3U  8U  f  1  ”  . 

C  c.  ~  2  i.  2.  £—•  ;yyr-  l  j-  1  op.  6  L. 

pc  3p .  3£.  M  rim  lm 

1=1  j=l  "l  ^i  m=l  J 


(A- 9 ) 


We  assume  that  deviations  in  p  ,  i  =  1,  2,  n,  are  not  correlated  with 

deviations  in  £ . ,  j  =  1,  2,  n.  Then  for  large  M  the  expression  in  brackets 

in  Eq.  (A-9)  is  taken  to  be  zero.  Hence  the  cross  term  is  zero. 

A  similar  argument  can  be  used  to  show  that  all  the  cross  terms  in  Eq.  (A-6) 
are  zero.  Thus  we  can  rewrite  this  equation  as 


T  +  T  +...+  T  +  Tf.  +  ...  , 

i/  p  £ 


P  v 


(A-10) 


where 


,3U  1  J 

W  M  \ 


(A— 11) 


T  =  (!^-)2  I  I  (6v  )2 

LV  '8vn '  M  .  v  0mj 
0  m=l 


(A-12) 
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T 


P 


I 

M 


M 

I 

ID- 1 


n 

I 

i  =  l 


au 

ap. 


6p 


ira 


i2 


(A- 13) 


and 


T 


i 


i 

M 


M 

I 

m=l 


au 


lm 


(A-14) 


The  mean  squared  deviations  in  Eqs .  (A-ll)  and  (A-12)  represent  variances. 
Thus  we  write 


T 


P 


3U  v 

ap 


and 


-  (  ajL  12 
"  1  9v„  J 


RvO 


(A-15) 


(A-16) 


where  a2n  is  the  variance  in  u„,  and  ct_2  is  the  variance  in  v„.  The  subscript 
K|JU  U  KVU  U 

R  indicates  that  the  launch  errors  are  taken  to  be  random. 


We  assume  that  the  deviation  6p.  from  the  mean  is  the  sum  of  a  "bias"  deviation 

lm 

6p_ .  and  a  random  deviation  Sp_.  .  The  6p_ .  are  correlated  from  zone  to 
Bim  rRim  Bim 

zone,  while  the  6p^m  are  uncorrelated.  Then  Eq.  (A— 13)  becomes 


T 


P 


1 

M 


M 

I 

m=l 


Z  3p .  (5pBim 

i=l  r  i 


6pRim) 


(A- 17) 
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Upon  expansion  of  the  square,  Eq.  (A-J7)  becomes 


M  n 


T  =  -  2 

p  M 


3U 


M  n 


2  §p7  6*W2  +  S  1  !  1 


m=l  i=l 


3U  .  12 

.  .  ,  3p.  ^PRim 

m=l  i=I  i 


(A- 18 ) 


where  the  omitted  cross  term  is  zero  because  the  6p_ .  are  not  correlated  with 

Rim 

the  6pBim' 

We  rewrite  Eq.  (A— 18)  as 


T  +  T 
lBp  Rp 


CA-19) 


where 


Bp 


M  n 


B  ^  1  2  fh  6pBim  l2 

m=l  i=l  ri 


(A-20) 


and 


Rp 


1  v  f  V  §U_  An  12 

M  2  1  2  3p.  6pRim  J 

m=l  i=l  l 


CA-21) 


After  expansion  and  rearrangement  of  sums,  Eq.  (A-21)  becomes 


V  =  ,2t  (ijr)2  f  r  \  * 


(A-22) 


where  the  omitted  cross  terms  again  are  zero.  In  Eq.  (A-22)  the  expression  in 
brackets  is  the  variance  computed  from  the  random  deviations  in  the  variable 

p^  appropriate  to  zone  top  i.  Thus 
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(A— 23 ) 


Rp 


2  (§»L)2 

l3p- 

1=1  i 


Rpi 


When  the  square  in  Eq.  (A-20)  is  expanded,  the  cross  terms  cannot  be  omitted 

because  of  the  correlation  of  the  6p„ .  from  zone  to  zone.  Equation  (A-20) 

Him 

becomes 


Bp 


n 

I 

i=l 


M 

I 

m=I 


<6pBim>2  I 


(A-24) 


n-1 

1 

i=l 


n 

I 

j  =  i+l 


au 

3p, 


au 


i 

M 


M 

1 

m=l 


6pBim  6pBjm  5 


If  is  the  bias  contribution  to  the  variance  in  p^ ,  then  Eq.  (A-24)  may 

be  written  as 


n 

1 

i=l 


l3p.j 


+  2 


n- 1 
I 

i=l 


n 

I 

j  =  i+l 


au 

aP, 


au 

apj  rij  aBpi  Bpj  ’ 


(A-25) 


where  r^j  is  the  correlation  coefficient  relating  errors  for  zones  i  and  j. 
We  assume  complete  correlation  for  the  bias  errors,  i.e.,  r^j  =  1.  Under  this 
condition  the  right  side  of  Eq .  (A-25)  reduces  to  the  square  of  the  sum  of 
terms  in  ogpi,  yielding 
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(A-26) 


! 

! 

| 

i 


TBp  -  <  £  f:  °Bpi  12 


From  Eqs .  (A-19),  (A-23),  and  (A-26),  we  have 


_  r  "  9U  ,o  ^  £  ,  dU  l2 

P  ~  1  °Bpi  1  jlj  1  3pt  °Rpi  1 


(A-27) 


A  similar  expression  can  be  obtained  for  : 


T{  *  1  sq  °Ki  '*  *  ii,  '  st-  °BCi 


3U 


,  ,  3U  , 

+  2.  [  zrzr-  dp,..  ) 


(A-28) 


where  °g^  and  are,  respectively,  the  bias  and  random  contributions  to 

the  variance  in  £  . 

From  Eqs.  (A-10) ,  (A- 15) ,  (A- 16),  (A-27),  and  (A-28),  we  express  the  variance 
in  the  East  component  of  ballistic  wind  as 


2  r  911  t~r  1 2  .  r  BD  ,2  . 

°U  =  [  9p^  RpO  ^  [  avQ  RvO  J 


*  1  S;  aBPi  ]2  +  ^  S:  aB£i  ]2  + 


♦  2  f  I ?-  ^  ]2  +  2  [  oD..  ]2  + 


9U 


3p .  Rpi 
i=l  l 


*  1  aT  RCi 

1=1  3i  ” 


(A-29) 
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APPENDIX  B.  PARTIAL  DERIVATIVES 


In  this  appendix  the  computational  forms  of  the  partial  derivatives  required 
by  the  error  models  are  given. 

With  reference  to  Figure  1,  the  following  relationships  can  be  written  for  the 
radiosonde  at  the  top  of  zone  i: 


x .  =  D .  sin  a . 

111 


(B-l) 


and 


y.  =  D.  cos  a. 
Ji  1  1 


(B-2) 


where  each  symbol  retains  its  definition  from  Section  2. 

Replacing  the  zone  index  i  with  the  launch  index  0  in  Eqs .  (B-l)  and  (B-2),  we 
can  evaluate  the  partial  derivatives  pertinent  to  the  launch  site.  For  example, 


sin  or„ 


(B-3) 


All  the  required  launch  partials  are  given  in  Table  XII. 

3Xj  8x, 

The  zonal  partial  derivatives  t— -  ,  -r —  ,  etc.,  i  =  1,  2,  . . . ,  n,  are  eval- 

9ti  9£i 

uated  assuming  a  spherical  earth.  To  achieve  this,  any  of  several  equivalent 
relationships  among  the  geometric  variables  may  be  used  as  a  starting  point. 
From  Reference  1  or  5,  we  write  the  distance  along  the  surface  of  the 
spherical  earth  as 
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R  COS  £ . 

D.  =  R  arccos  (  -=-  - - -  )  -  Re.  ,  (B-4) 

1  K  +  z .  1 

1 

where  R  is  the  mean  radius  of  the  earth  and  z^  is  the  altitude  at  the  top  of 
zone  i.  The  elevation  angle  e is  expressed  in  radians. 

From  Eqs.  (B-l),  (B-2),  and  (B-4) ,  we  have 


R  cos  £  . 

x.  =  [  R  arccos  (  — = - -  )  -  Re.  ]  sin  a.  ,  (B-5) 

i  x  K  +  z  i  i 

i 

and 

R  cos  £ . 

y.  =  [  R  arccos  (  g- - -  )  -  Re.  ]  cos  a.  ,  (B-6) 

1  K  t  z ,  1  1 

1 

The  required  partial  derivatives  of  and  y^  with  respect  to  z^,  e^,  and 
are  obtained  from  direct  dif ferentation  of  Eqs.  (B-5)  and  (B-6).  The  results 
are  given  in  Table  a"'I.  In  order  to  display  the  results  compactly,  it  is 
convenient  to  define  the  following. 


«i  = 


R  +  z. 


(B— 7  ) 


and 


d».  =  arccos  (Q.  cos  e.) 

r  i 


(B-8) 
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3z .  9z . 

The  partial  derivatives  gjp  and  g^—  are  required  by  the  RADAR  model.  From 

i  i 

Reference  6  the  slant  range  is  given  by 


S. 

l 


[(z.  +  R)2  +  R2  cos2  e.]1'* 


R  sin  £ . 

i 


(B— 9 ) 


We  could  solve  this  equation  for  z^  and  proceed  to  evaluate  the  needed  partials. 
However,  it  is  somewhat  easier  to  pursue  the  procedure  outlined  below. 

We  define  the  function  f(S.,  £.,  z.)  to  be 

l  ’  i  l 


f  (S .  ,  £.,  z.)  =  l(z.  +  R)2  -  R2  cos2  e.]"4 

i  ’  i  ’  l  l  i 


R  sin  £ ^  (B-10) 


Then  from  Eqs.  (B-9)  and  (B-10), 

f  (S  . ,  £  .  ,  z  .  )  =  0 

l  l  i 


(B- 11) 


The  altitude  z^  is  taken  to  be  dependent  on  and  £^. 

Given  the  condition  stated  in  Eq.  (B— 11),  a  theorem  of  partial  differentiation 
allows  us  to  write 


3z. 

l 

as. 


3f  .  3f_ 

as .  '  dz. 

i  i 


(B- 12) 
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when  is  not  equal  to  zero.  See  Reference  7,  Chapter  5,  for  example,  for 

i 

a  discussion  of  this  theorem. 


The  needed  partial  derivatives  of  f  are  evaluated  using  Eq.  (B-10).  Then 
3z .  3z . 

gg-  and  are  obtained  from  Eqs .  (B-12)  and  (B-13),  respectively.  The 

i  i 

results  are  displayed  in  Table  XII. 


Table  XII.  Partial  Derivatives 


3D„ 


=  sm  a. 


3yc 

3D„ 


=  cos  a. 


3xc 

3a, 


D0  cos  aQ 


3yC 

3a, 


Do  sin  ao 


3x .  Q. 2  cos  s .  sin  a . 
i  _  i  i  i 


dz . 
i 


sin  0 . 

i 


dy .  Q.2  cos  £.  cos  a. 

i  i  i  i 


dz  . 

i 


sin  0 . 


3x . 
i 

de. 

i 


Q  .  sin  6 . 

R  sin  a.  [1  -  — ^ — 7 — r — - 

i  sm  0. 


dji 

3e . 
i 


-  R  cos  a^  [ 1 


3x . 
i 

3a. 

i 


R  (()) .  -  e  . )  cos  a . 
i  i  i 


3yi 

3a. 


R  (0 .  -  £  .  )  sin 
i  l 


3z . 
_ i 

3S. 

i 


sin  0^ 


3z  . 

i 

3  e. 

i 


R  sin  (0^  -  £ ^ ) 


where  R  =  mean  radius  of  the  earth 


Q- 


R 

R  +  z. 
i 


0.  =  arccos  (Q.  cos  £.) 

i  i 
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